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Abstract. - We study spin-wave interactions in quantum antiferromagnets by expressing the 
usual magnon annihilation and creation operators in terms of Hermitian field operators rep- 
resenting transverse staggered and ferromagnetic spin fluctuations. In this parameterization, 
which was anticipated by Anderson in 1952, the two-body interaction vertex between staggered 
spin fluctuations vanishes at long wavelengths. We derive a new effective action for the stag- 
gered fluctuations only by tracing out the ferromagnetic fluctuations. To one loop order, the 
renormalization group flow agrees with the nonlinear-(T-model approach. 



The spin-wave approach to quantum Heisenberg antiferromagnets, which was pioneered 
by Anderson [1] and Kubo [2] more than half a century ago, is still one of the most powerful 
methods to calculate the low-temperature properties of ordered magnets. In this approach, the 
components of the quantum-mechanical spin-operators Si are expressed in terms of canonical 
boson operators bi using the Holstein-Primakoff [3] or the Dyson-Maleev [4] transformation. 
The spin-S' antiferromagnetic Heisenberg Hamiltonian 



H 



r, ^ JijSi ■ Sj , (1) 



can then be written as a bosonic many-body Hamiltonian H = —DNJS^ + H2 + -ffint, where 
the quadratic part is 

H2^SJ2 J^o + b]b, + 6.6, + 6l6]] . (2) 

The spins are assumed here to be localized at the sites of a ZJ-dimcnsional hypercubic 
lattice with lattice spacing a. The lattice is bipartite, with sublattices labelled A and B. We 
also assume = J > for all pairs of nearest neighbors and Jtj = otherwise. 

The interaction J?int involves at least four boson operators and higher powers of the small 
parameter l/S", so that for large S it is reasonable to treat -ffint perturbatively. However, 
as we discuss in more detail below, in momentum space the vertices of -ffint have a compli- 
cated non-analytic structure for small momenta [5-7] , so that the expected supression of the 
interaction between long-wavelength Goldstone modes [8] is not manifest in this approach. 
On the other hand, in the calculation of physical quantities that are dominated by staggered 
spin-fluctuations the leading momentum-dependence of the vertices eventually cancels [9]. 
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Note that the l/S'-expansion for quantum ferromagnets is not plagued by this problem: in 
this case the two-body interaction between ferromagnetic magnons vanishes quadratically for 
small momenta, reflecting the supression of the effective interaction between long-wavelength 
Goldstone-modes [fO]. 

Alternatively, we may work with the nonlinear sigma model (NLSM), which is believed to 
describe the long- wavelength and low-energy physics of quantum Heisenberg antiferromagnets 
in the so-called renormalized classical regime [11] (for a careful discussion and a derivation of 
an effective model for the short wavelength regime see [12]). In imaginary time r and at finite 
temperature T = 1/ (3 the action of the NLSM is 



?nlsm[J^] = Y J drld^r 



{d,ny + c,'{drny , (3) 



where the unit vector r2(r,r) represents the slowly fluctuating staggered magnetization, po 
and Co are the spin stifl^ness and the spin-wave velocity at T = 0. Assuming local staggered 
order in z-direction, we may resolve the constraint ft^ — 1 by writing the z-component of ft 
as i/l — « 1 — in^ -I- . . ., where the two-component vector 11 represents the transverse 
staggered spin fluctuations. The interaction vertices generated in Eq. (PJ by expanding the 
square root involve two derivatives, so that their Fourier transform vanishes quadratically for 
small wave- vectors or frequencies, in contrast to the vertices in the usual l/S'-expansion [5,6]. 

Whether it is possible to parameterize the 1 / S-expansion such that the weakness of the 
interaction between long-wavelength staggered spin fluctuations is manifest is a long-standing 
unsolved problem in spin-wave theory, which we shall solve in this work. Let us therefore 
recall the usual diagonalization procedure of the quadratic spin- wave Hamiltonian ||2Jl. First 
of all, we introduce the Fourier components of the operators bi in the sublattice basis, defining 
bi = (2/A^)i/2 J2k e^^-^'^Ak if n belongs to sublattice A, and 6; = {2/Nfl'^ Y.k e'^ ^'^Bk if 
belongs to sublattice B. Here and below the fc-sums are over the reduced Brillouin zone. The 
complete diagonalization of H2 is then achieved with the help of a Bogoliubov transformation, 

.Blk)-[ j Ul J ' 

where Uk = + l]i/2 and Vk = 2-i/2[e^i - l]i/2 ^jth = [1 - Here 

7fc = D^^ cos(fc^a), where ~ e^-k are the components of k in the direction of the unit 
vectors e^, fx— 1, . . . ,D. Then we obtain H2 = —NDJS + H2, where 

k 

with the magnon dispersion Ek = 2DJSek. The Bogoliubov transformation is not unique: 
any rotation that mixes the operators ak and (3k will also diagonalize H2. To make contact 
with the H-field in the NLSM, we choose ^fcc, — 2^^/'^{ak + crf3k), where a — ± labels 
the symmetric and antisymmetric combinations. We then express the magnon annihilation 
operators '^ka in terms of two Hermitian field operators Uka and ^ka as follows, 

^-fca = PaiXof^VEky/^EkUka+iXo'^ka] , (6) 

where xo — {4:DJa^)^^ = Pq/cq is the classical uniform transverse susceptibility. The phase 
factors = —i and p_ — 1 are introduced for later convenience and V = N is the volume. 
One easily verifies that [Hfeo-, ^fcv] = 'iVSk,-k'Sa-,(T' , so that Uka and ^ka are canonically 
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conjugate bosonic field operators, with likal'^ka corresponding to the position/momentum 
operators. These operators are simply related to the Fourier components of the total spin 
Sk = 7V-i/2^,e-^'= '--5, and the staggered spin S.t.k = N-^/^ J2^e-'''-^%,S,, where Q = 1 
for Ti £ A and (^j = — 1 for e i?. For large S we find to leading order 



and 



Stt,k « N-^'\S/a^)\ktVk+, 
Slk « N-'/'{S/a^)X^Uk^, 



(7a) 
(7b) 

(8a) 
(8b) 



Here Afc = [uk + Wfc](efc/2)^/^ = [(1 + 7fc)/2]^/^ approaches unity for fc ^ 0. Hence, our 
Hermitian field operators can be identified physically with the suitably normalized transverse 
components of the staggered and total (ferromagnetic) spin for large S. Writing $ — (<i>-|_, <&_) 
and n = (n+, n_) and using Eq. we find that H2 takes the form 

^2 - y}_,[ +^n_..n.j . (9) 

The parameterization of the spin fluctuations introduced by Anderson in his early work on 
the antiferromagnetic ground state [1] differs from our Eq. ^ only in the choice of normal- 
ization factors. Subsequently, this method has been far less popular than the approach based 
on the usual bosonic annihilation and creation operators and has been largely forgotten. Nev- 
ertheless, the Hermition operator approach has many advantages: (a) the suppression of the 
effective interaction between antiferromagnetic magnons at long wavelengths is manifest; (b) 
the relation between the Holstein-Primakoff bosons and the continuum fields 11 in the NLSM 
can be made precise; (c) a new effective action for the transverse staggered spin fluctuations 
can be obtained by eliminating the ferromagnetic fluctuations in a path integral formulation. 
The resulting interaction vertices differ from those of the NLSM, because the l/S-approach 
contains scattering processes which are not related to the constant length constraint of the 
spin vector. At one-loop order it turns out that the two-body interaction does not renormal- 
ize the Gaussian part of the effective action, so that the quantum critical point separating 
the renormalized classical from the quantum disordered phase can be directly related to the 
anomalous dimension of the Il-field; (d) the Hermitian operator approach is most convenient 
to discuss spin-waves in finite-size antiferromagnets. 

We now explain the above points in some detail, beginning with (d). As first pointed 
out by Anderson [1], for finite N the fc = term in Eq. should be treated with special 
care. Defining the dimcnsionless operators P = {SN)^^^^^o and X = {S/Ny^^a~-^Ilo 
whose components satisfy [Xo-,^^'] = iSa^a', and using the fact that for periodic boundary 
conditions Ek=o = 0, the contribution from the fc = term to Eq. ^ is H2 = ^ with 
m = ADJS. Obviously, the spectrum of H2 is continuous; the ground state is the zero- 
momentum state satisfying (P^) = 0. However, for any state with {X) = (P) = we have 
(P^)(X^) > 1 by the uncertainty principle, so that {X'^) = 00 in the ground state. The 
staggered magnetization Mst ~ Ci'S'f ) can be written as 



M,t = N{S + 1/2) - (P2) - (X2) - — Y,{h^-k ■ + ./fc 'n-fe • tlk) , (10) 

k=iO 
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where ff, — a^/{SX^). In the ground state of H2 the zero mode gives rise to an infinite 
correction to Mgt- At first sight, it appears that our approach is inconsistent. The reason why 
the spin-wave approach can still be useful in finite systems has been discussed by Anderson 
[1]: suppose we prepare a finite-size antiferromagnet by some external field in a minimal 
uncertainty wave-packet with (P^) = {X'^) = 1. The zero-mode contribution to Mst is 
then of relative order 1/N, and can be ignored for large N. The validity of the spin- wave 
approach becomes then a dynamical problem: from quantum mechanics we know that under 
the time evolution governed by H2 — ^ the width of the Gaussian wave-packet increases 
as {X^)t — {X^)o + ^{P'^)o so that the spin-wave approximation breaks down at time t « 
\fNm = / (4:0,78), which for macroscopic systems exceeds the time scale of experiments. 

We now go beyond Ref. [1] and consider interactions between spin-waves. Within the 
Holstein-Primakoff transformation, the leading interaction correction is given by the following 
bosonic two-body Hamiltonian 

H4 = X! "^"-i [^"^i^j + iT'ibibj + biUjbj + b\b^.nj + fofr^ife]] , (U) 

where rii = h\hi. It is now straightforward (although quite tedious) to express in terms 
of our Hermitian field operators. Since we shall later use the phase space path integral [13] 
to eliminate the ferromagnetic fiuctuations, we symmetrize H/^ whenever the vertices involve 
non-commuting operators [14]. The final result is 

" + § + ^ 4i+fc.+fc3+fc4,o[r™(fei,fc2,fc3,fc4)rifei • rife.nfe, • 

fci ,...,A;4 

+ /orf "(fci, fc2, fea, fc4) \{^k,Ak,a.%,af^k,.] 
a 

+fSrf{ki,k2,k3, ki) ^ $fc,a^fc..nfc3,-^nfe,,_, 
+/2r*n(fci,fc2,fe3,fc4){<i'fei+,nfc.+}{nfc3_,l>fc,.} 

+ /o4r**(fci,fc2,fc3,fc4)*fei -^fc.^fca •*fe4] . (12) 

where {A, _B} = AB + BA denotes the anticommutator, = —{DJ/A)[N + X]fc7fe]' ^"^^ 
fo = fk=o — / S. Writing 71 = 7^^ etc., the properly symmetrized vertices in Eq. H12|) are 

r™(fci,/c2,fc3,fc4) = (l/8)AiA2A3A4[7i+72 + 73 + 74 - 2(7i+2+73+4)] , (13a) 
rfn(fci,/c2,fc3,fc4) = (l/4)(AiA2)-^A3A4 [-71 -72+73+74-2(71+2+73+4)] ,(13b) 
Tf'{ki,k2,k3, ki) = (l/4)(AiA2)-^A3A4 [3(-7i - 72 + 73 + 74) - 2(71+2 + 73+4) 

+4(71+3 + 72+4)] , (13c) 
Vl^'ik^MM. ki) - (l/8)Ai IA2A3A4 1 [71 - 72 - 73 + 74 - 4(71+3 + 72+4)] , (13d) 
r**(fci,fc2,fe3,fe4) = (l/8)(AiA2A3A4)-^[-7i-72 -73 -74 - 2(7i+2+73+4)] .(13e) 

The important point is now that the vertex r'^'^ associated with the direct interaction between 
staggered spin fluctuations vanishes quadratically for small momenta, 

T^^{ki,k2,k3, fc4) = (aVl6-D)[fci +ki + kl + kl+ 4fci • ^2 + 4fc3 • ^4] + O(fcf) , (14) 
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while the other vertices approach finite limits, 

r|n(o) = r,*n(0) = r**(0) = -l , and r*n(o) = 1 . (15) 

Hence, in our parameterization the bare interaction between the staggered spin fluctuations 
is manifestly suppressed at long wavelengths, as seen from Eq. H14|l . General symmetry argu- 
ments [8] suggest that this suppression survives when we take into account the renormalization 
of the staggered spin fluctuations by the ferromagnetic ones. We have explicitly verified this 
within the Hartree-Fock approximation, where the two-body part in Eq. H12(l is replaced by 
a one-body Hamiltonian. Using {fl^ka'^ka) = V^(2xo£'fe)~\ {^-ka^ka) = VxoEk/l, and 
(Jl-ka'^ka) = *^/2, we find that the Hartree-Fock approximation amounts to the replace- 
ment E'^ + jgH2, where E'^ differs from i?4 given above by a constant of the order of 
unity, and C = 1 — {2/N) e^. It follows that the leading l/S'-correction to the magnon 
dispersion can be taken into account in the free spin-wave Hamiltonian H'2 given in Eq. Q 
by replacing xo ■^xXo a-nd Ek — > ZcEk, with the renormalization factors = [1 + ^]~^ 
and Zc = 1 + ^ . Note that this simply leads to an overall rescaling of the magnon dispersion, 
but does not change its wave- vector dependence [15]. 

In contrast, the vertices encountered in the usual spin- wave theory, based on the magnon 
operators ak and /3fc defined in Eq. Q , have a more complicated structure. Using the Dyson- 
Maleev transformation (the vertices derived from the Holstein-Primakoff transformation are 
linear combinations of the Dyson-Maleev vertices [5]), the interaction part takes the form 
[16,17] 

^DM ^ ^ ^ {yW4(/3t/?|/33/34 + a^4«ia2)-2y/|^4(at/34aia2 + 4/314/^3) 

ki ,...,^4 

+ '^^1234("3"4/3l/32 + Olia2P-iP^5ki+k2M3+ki ■ (16) 

For simplicty, we have introduced the notation ai = aki etc. In the limit that all wave vectors 
are small, the Dyson-Maleev vertices v/234 = ^^■'H^i: ^2, ^3, ^4) behave as [16] 

yU) _ 1 / |fcil|fc2| / fci ■ fc2 \ ^ 

with ^1 = ^2 = ^5 = 1 a^^d ^3 = ^4 = —1. Obviously, these vertices do not vanish for 
small momenta. Moreover, the long-wavelength limits are direction-dependent and the non- 
analytic prefactor can potentially give rise to divergences in perturbation theory. Although 
the weakness of the underlying spin- wave interaction is not aparent, it was shown in [6] 
that to order 1/5^ the divergences cancel in a l/S'-expansion as a consequence of total spin 
conservation [9]. Comparing Eq. ifTTjl to the equivalent Eqs. H14I15|I of the Hermitian operator 
formulation, the much simpler structure of the latter becomes clear. This simpler structure is 
a direct consequence of using physically transparent operators. 

In the conventional approach based on the operators ak and Pk, it is very cumbersome 
to calculate higher order terms in a l/S'-expansion [5,6]. Instead of performing a complete 
l/S'-expansion, we shall here only concentrate on the long wavelength antiferromagnetic fluc- 
tuations. Within the Hermitian field parameterization these are clearly separated from the 
ferromagnetic ones. Moreover, the weakness of the interaction between long wavelength stag- 
gered fluctuations is explicit in this approach. We thus simply eliminate the ferromagnetic 
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degrees of freedom represented by the ^-operators and work directly with the effective ac- 
tion S'cfr[n] of the staggered fluctuations. Formally, S'off[n] can be defined as a phase space 
path-integral [13,14], g-'^-ffM = J 'D[^]e-^^^'^\ where the Euclidean action S'[n, *] is a 
functional of quantum fields nfc(r) and $fe(T) associated with the field operators. Within the 
Gaussian approximation we obtain 

S[U, *] » i / [xo • + XoElU^K ■ Uk - ujni^-K ■ Uk - n K ■ ^k)] , (18) 

where K = (k, iujn) is a collective label for wave- vector and bosonic Matsubara frequency w„, 

the Fourier transformed fields are Hk = Jo dre^'^^'^IlkiT), and — {f3V)~^ J^ku^- The 
corresponding Gaussian effective action for the staggered fluctuations is 

S,s[n] ^^J^{El+ ujDu^k ■ Uk . (19) 

In the long- wavelength limit ~ co|fe|, so that Eq. (|19|l has precisely the same form as the 
Gaussian approximation of the NLSM in Eq. (PJ . Because we have explicitly constructed the 
H-field in Eq. (|19|l from the Holstein-Primakoff bosons bi , we have succeeded to give a mathe- 
matically precise relation between the two different parameterization of the spin fluctuations. 

The leading interaction correction to iSeff [H] is due to the vertex in Eq. (|12|l involving 
four Il-operators. Retaining only this term we obtain in the long-wavelength limit with the 
help of Eq. ifT^ at T = 0. 

5cff[n] « J-JdxoJ d"x[{d^ii)^ + idoUf + l[2{u-d^uf-u^d^u)^]] , (20) 

where x = Apr and xq = AqCqt are dimensionless space-time variables, go = Aq~^cq/po is the 
usual [11] dimensionless coupling constant of the NLSM, and Aq ~ 1/a is an ultraviolet cutoff. 
The quartic interaction in Eq. H2U|I does not agree with the quartic term in the perturbative 
expansion of the action of the NLSM [11]. In the latter case, the term in the square braces 
of Eq. (|20|l is replaced by 4(11 • 9/^11)^ and there is an additional counter-term due to the 
expansion of the functional 5- function enforcing the constraint fl^ — 1. The reason for this 
difference is that in the NLSM all interactions arise from the non- linear constraint ft^ — 1. In 
contrast, the l/S'-approach contains dynamic interactions neglected in the NLSM, but ignores 
the constraints associated with the finiteness of the spin Hilbert space which according to 
Dyson [4] amount to irrelevant kinematic interactions. 

It turns out that within a one-loop momentum shell renormalization group (RG) approach 
[10,11] the quartic interaction in Eq. H2()|l does not renormalize the Gaussian part H19|l of our 
effective action at all. A similar cancellation is known to happen also in the corresponding 
ferromagnetic Heisenberg model [10]. If we apply the usual momentum scale RG procedure 
consisting of mode elimination and rescaling [10, 11], then the dependence of the running 
coupling gi on the logarithmic RG flow parameter I is at one loop entirely due to the rescaling 
step. The resulting one-loop RG flow equation is simply 

dm^-iD-i-vDgi- (21) 

Here rj^ is the running anomalous dimension of the H-field (not to be confused with the 
usual critical exponent ij), which is related to the interaction-dependent part of the field 

1 /2 

rescaling factor as usual, 77^ — —dilnZi. The simplest way to calculate ryf is from the 
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transformation properties of the staggered magnetization il/gt given in Eq. H1U|1 under the 
momentum scale RG [10]. By demanding that the scahng of the fluctuation correction in 
Eq. H10|l is consistent with the leading term N{S + 1/2) (which scales as the volume), it is 
easy to show that -qf = ^f-gi, where Kjj is the surface of the unit sphere in D dimensions 
divided by (27r)-°. The resulting RG equation for gi, 

digi^{l~D)gi + ]^KDgh (22) 

is of course identical to the equation derived in Ref. [11]. Interestingly, in our formulation 
the quantum critical point separating the renormalized classical from the quantum disordered 
regime can be characterized by lim/_,oo vij ^ D — 1. 

In summary, we have developed a new parameterization of the spin- wave expansion which 
is clearly superior to the conventional formulation based on the magnon operators au and 
f3k in Eq. Q). In our parameterization the physical meaning of the degrees of freedom is 
transparent, the weak interaction between long- wavelength staggered magnons is manifest, 
the subtleties of the spin-wave approach in finite-size systems can be easily discussed, and 
the effective theory for the staggered fluctuations can be derived. With our method it is also 
straightforward to discuss the effect of interactions between the zero modes on the flnite-size 
spin- wave spectrum, and to study spin- wave interactions in more complicated models involving 
various anisotropics or external magnetic fields [18]. 

* * * 
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